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1.0 SUMMARY 


Program L219 (EQMOD) is structured as eight overlays, one main, four primary, and 
three secondary overlays. Input into the program is made via cards and magnetic files 
(tapes or disks). Output from the program consists of printed results and magnetic files 
containing modified copies of the input matrices. 

Although L219 (EQMOD) serves as a module of the DYLOFLEX system, it can be 
operated as a standalone program. Subroutines used by L219 include routines embedded 
in the program code and routines obtained from the standard FORTRAN, DYLOFLEX, 
and FLEXSTAB libraries. 


2.0 INTRODUCTION 


Program L219 (EQMOD) can be used as either a standalone program or as a module of a 
program system called DYLOFLEX (see fig. 1) which was developed for NASA under 
contract NASl-13918 (ref. 1). Because of the DYLOFLEX contract requirements 
developed in reference 2, a program was needed to modify the equations of motion 
matrices generated by L217 (EOM) of reference- 3 and the load equation matrices 
generated by L218 (LOADS) described in reference 4. The matrices are modified to 
include: 

• Scalar multipliers 

• Replace or increment individual matrix elements 

• Add sensor equations to the equations of motion 

• Add the definition of the active control system to the equations of motion 

• Replace the rigid body stability derivatives in the equations of motion with those 
calculated by FLEXSTAB (ref. 5) or other external means 

• Transform the equations of motion and load equations from the inertia axis system 
to the body axis system 

• Prepare the matrix coefficients in a form usable by the Linear System Analysis 
program QR (ref. 6). including: 

• Equations of motion with and without Wagner lift growth functions 

• Equations of motion and load equations combined for a time history solution 




Figure 1.- DYLOFLEX Flow Chart 









The objective of this volume is to aid those persons who will maintain and/or modify the 
program in the future. To meet this objective, the following items are defined in detail: 

• Program design and structure 

• Overlay purpose and. description 

• Input, output and internal data base descriptions 

• Test cases 

The program was designed, coded, and tested according to the DYLOFLEX 
programming specifications . 1 


! R. E. Clemmons: Programming Specifications for Modules of the. Dynamic Loads Analysis System 
to Interface with FLEXSTAB. NASA contract NAS1-1391B. BCS-G0701 (internal document) , 
September 1975. 
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3.0 PROGRAM DESIGN AND STRUCTURE 


Program L219 (EQMOD) has been constructed as an overlay system consisting of a 
main overlay, four primary and three secondary overlays (see fig. 2). 


Main overlay (L2 19,0,0) L219vc 

Primary overlay (L219,l,0) RDCRDS 

Secondary overlay (L219,l,l) RDEOM 

Secondary overlay (L219,l,2) RDLOD 

Secondary overlay (L219,l,3) RDQR 

Primary overlay (L219,2,0) EOMMOD 

Primary overlay (L219,3,0) LODMOD 

Primary overlay (L2 19,4,0) QRMOD 


The input and output of EQMOD are displayed in figure 3. The file scratch random file 
SCRAND is used for temporary storage of data by EQMOD. The other files 
communicate with programs outside of EQMOD. 

The main overlay L219vc initializes the program’s default values and the scratch 
random file, determines which primary overlays are to be executed, and aids 
communication between the overlays through labeled common blocks. The characters 
v and c in the program name stand for version and correction respectively (see sec. 3.1). 

The 1,0 primary overlay RDCRDS reads and edits common input data, writes 
diagnostics when errors are detected, and determines which secondary overlays are to 
be executed. In addition, if FLEXSTAB rigid body stability derivatives are requested. 
RDCRDS will read the stability derivatives from file SDSSTP. The data read from cards 
and SDSSTP are stored in labeled common blocks and in arrays written onto the scratch 
random file SCRAND. RDCRDS calls RDEOM, RDLOD, and RDQR overlays to perform 
portions of its tasks. 

The 1,1 secondary overlay RDEOM reads and edits all Equations of Motion (L217) input 
data, writes error diagnostics when an error occurs, and stores edited data on the 
scratch random file SCRAND. 

The 1,2 secondary overlay RDLOD reads and edits all Load Equations (L218) input 
data, writes error diagnostics when an error occurs, and stores edited data on the 
scratch random file SCRAND. 

The 1,3 secondary overlay RDQR reads and edits all Linear Systems Analysis QR 
program input data, writes error diagnostics when an error occurs, and stores edited 
data on the scratch random file SCRAND. 
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Figure 2. — Overlay Structure of L219 (EQMOD) 











EQEOM EQLOD QRTAP 


The input and output magnetic files have 
variable user specified names. 


Figure 3.— Input/output of L219 (EQMOD) 
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The 2,0 primary overlay EOMMOD modifies all of the matrix coefficients from the 
Equations of Motion (L217) program according to the instructions specified on input 
cards. 

The 3,0 primary overlay LODMOD modifies all of the matrix coefficients from the Load 
Equations (L218) program according to the instructions specified on input cards. 

The 4,0 primary overlay, QRMOD generates rooting and time solution matrices for the 
Linear Systems Analysis, QR, program. 

Each overlay is discussed in detail in succeeding sections. Included for each overlay are: 

1. The overlay’s purpose. 

2. The overlay’s analytical steps. 

3. The input/output devices used. 

4. A macro flow chart. 

5. Table of subroutines called. (Note: all subroutines have only one entry point.) 

3.1 MAIN OVERLAY (L219,0,0)-L219vc 

The main overlay of L219 (EQMOD) is L219vc, where v is a letter indicating the 
program version, and c is an integer number indicating the correction number which 
applies to the v version. 


Purpose of L219vc 

L219vc performs certain bookkeeping tasks, directs the execution of the primary 
overlays, and aids communication between primary overlays via labeled common blocks. 

The macro flow chart of this overlay is shown in figure 4. The subroutines called are 
displayed in table 1. 


I/O Devices of L219 

L219 reads a data card (card set 1.0) and writes diagnostics on the output file if errors 
are encountered. All other I/O accomplished by L219 (EQMOD) is done in lower level 
overlays. 


3.2 PRIMARY OVERLAY (L219,1,0)-RDCRDS 
Purpose of RDCRDS 

The primary overlay RDCRDS is called to read and edit common input data, print a 
diagnostic when an error is detected, and call secondary overlays to read L217 (EOM), 
L218 (LOADS), and QR input data. 


7 



Figure 4. -Macro Flow Chart of Overlay (L219, 0, 0)-L219vc 
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Table 1.— Routines Called by L219vc 


OVERLAY (L2 19,0,0) 
PROGRAM L219vc 


RDCRDS (Overlay L219,1,0) 
EOMMOD (Overlay L219,2,0) 
LODMOD (Overlay L219,3,0) 
QRMOD (Overlay L2 19,4,0) 
CLOSMS * 

DATE * 

FETAD + 

FETDEL ♦ 

IPQL + 

KEYWRD calls EOF * 
OPENMS * 

OVERLAY * 

PRGBEG + 

PRGEND ♦ 

RETURNF + 


+ indicates DYLOFLEX library routine 

* indicates 6600 system library routine 


I II II II I II 


Analytical Steps of RDCRDS 

RDCRDS reads a card, extracts the keyword code, and performs one of the following 
tasks before reading another keyword card. 

• Store title card 

• Define problem size 

• Define output tape name and matrix positions 

• Define print option data 

• Set symmetric/antisymmetric indicator 

• Call program RDEOM to read and edit EOM data 

• Call program RDLOD to read and edit LOADS data 

• Call program RDQR to read and edit QR data 

• Print error diagnostics when errors are detected 

The macro flow chart of this overlay is shown in Figure 5. The subroutines called are 
displayed in table 2. 


I/O Devices of RDCRDS 

RDCRDS reads card sets 1.0 through 19.0 and echo prints the input data and 
diagnostics for all errors detected. 

3.2.1 SECONDARY OVERLAY (L219,l,l)-RDEOM 

Purpose of RDEOM 

The secondary overlay RDEOM is called to read and edit all Equations of Motion L217 
(EOM) input data, write error diagnostics when an error is detected, and store edited 
data on scratch random file SCRAND for the EOM equation matrix modifier program. 

Analytical Steps of RDEOM 

RDEOM has the following three steps: 

1. Decode $EOM card for input tape name, matrix file position, and null matrix 
indicators. 

2. Read a card, find keyword code, and perform one of the following tasks: 

• Call RDDERS or RDDERA to read, edit, and calculate derivatives for 
symmetric or antisymmetric analysis, respectively 

• Call RDSEN to read and edit all sensor data. 

• Call RDSAS to read and edit all active control system definition data 

• Call RDSCAL to read and edit all matrix scalar data 
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Table 2.— Routines Called by RDCRDS 


OVERLAY (L219, 1,0) 
PROGRAM RDCRDS 


RDEOM (Overlay L219,1,1) 
RDLOD (Overlay L219,1,2) 
RDQR (Overlay L219,1,3) 
KEYWRD 
NAMFIL + 

OVERLAY * 


+ indicates DYLOFLEX library routine 
* indicates 6600 system library routine 
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• Call RDINRE to read and edit all matrix element replacement and increment 
data 

• Read and edit body axis symmetric or anti-symmetric data 

3. Print diagnostic messages for all errors detected. The macro flow chart for RDEOM 
is in figure 6. The subroutines called are displayed in table 3. 

I/O Devices of RDEOM 

RDEOM reads card sets 7.0 through 13.0. It prints the calculated values with 
descriptive captions. 

RDEOM reads the header matrix from the file EOMTAP if the file was generated by the 
DYLOFLEX program L217 (EOM) (ref. 4). 

RDEOM writes on the scratch random file SCRAND all the edited EOM input data, the 
edited derivatives, and the active control system definition data. 
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ENTER 


Decode $EGM card 


EOMTAP 


Read header matrix 
from LODTAP 
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DERIVATIVE, 
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edit, and calculate derivatives 
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RETURN 


Figure 6. —Macro Flow Chart of Overlay ( L219 , 1 , D—RDEOM 
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Table 3.— Routines Called by RDEOM 


OVERLAY (1.219,1,1) 
PROGRAM RDEOM 


RDEOM calls FETAD + 

FETDEL + 

KEYWRD 
NAMFIL + 

RDDERA calls CLRTAB $ 

DMPTAB $ 
FETAD + 
FETCHM + 
FETDEL + 
INPTAB $ 
PRNTM 
WRITMS * 

RDDERS calls CLRTAB $ 

DMPTAB $ 
FETAD + 
FETCHM + 
FETDEL + 
INPTAB $ 
PRNTM 
WRITMS * 
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Table 3.— (Concluded) 


RDINRE calls IRQL + 

KEYWRD 

RDSAS calls KEYWRD 

WFITMS * 

RDSCAL calls KEYWRD 


RDSEN calls KEYWRD 

NAMFIL + 
WRITMS * 


READTP + 


WRITMS * 


$ indicates a FLEXSTAB library routine 
■f indicates DYLOFLEX library routine 
* indicates 6600 system library routine 
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3.2.2 SECONDARY OVERLAY (L219,l,2)-RDLOD 


Purpose of RDLOD 

The secondary overlay RDLOD is called to read and edit all Load Equations L218 
(LOADS) input data, write an error diagnostic when an error is detected, and store 
edited data on scratch random file SCRAND for the LOADS equation matrix modifier 
program. 


Analytical Steps of RDLOD 

RDLOD has the following four steps: 

1. Check if maximum number of LOADS sets processed 

2. Decode $LOADS card for input tape name, matrix file position, number of output 
loads, and null matrix indicators 

3. Read a card, find keyword code, and perform one of the following tasks: 

• Call RDSCAL to read and edit all matrix scalar data 

• Call RDINRE to read and edit all matrix element replacement and increment 
data 

4. Print error diagnostic when an error is detected 

Figure 7 is a macro flow chart of RDLOD. Table 4 displays the routines called by 
RDLOD. 


I/O Devices of RDLOD 

RDLOD reads card sets 15.0 through 17.0. 

If the loads input magnetic file, LODTAP was generated by the DYLOFLEX program 
L218 (LOADS) (ref. 4), the header array is read by RDLOD to determine which load 
matrices are available and their size. 

The input data and diagnostics are printed with descriptive labels. The edited loads 
input data is saved on the scratch random file SCRAND. 

3.2.3 SECONDARY OVERLAY (L219,1,3)-RDQR 

Purpose of RDQR 

The secondary overlay RDQR is called to read and edit all Linear Systems Analysis 
(QR) input data, write error diagnostics when an error is detected, and store edited data 
on the scratch random file SCRAND for the QR matrix generation overlay, QRMOD. 
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Table 4.— Routines Called by RDLOD. 


OVERLAY (L2 19,1,2) 
PROGRAM RDLOD 


RDLOD 

calls 

FETAD + 






FETDEL 

+ 





KEYWRD 






NAMFIL 

+ 





RDINRE 


calls 

IRQL + 






KEYWRD 



RDSCAL 


calls 

KEYWRD 



READTP 

+ 





WRITMS 

* 




+ indicates DYLOFLEX library routine 

* indicates 6600 system library routine 
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Analytical Steps of RDQR 


RDQR program has the following four steps: 

1. Check if the maximum number of QR sets has been processed 

2. Decode $QR card for output tape name and matrix positions 

3. Read a card, find keyword code, and perform one of the following tasks: 

• Read the EOM input tape name, matrix position, and null matrix indicators 

• Set Wagner function indicator and read Wagner functions 

• Set rooting function indicator 

• Set time solution indicator 

• Read the LOADS input file name, matrix position, number of loads, and null 
matrix indicators 

4. Print error diagnostics when errors are detected. 

Figure 8 is the macro flow chart of RDQR. Table 5 displays the subroutines called by 
RDQR. 

I/O Devices of RDQR 

RDQR reads cards 18.2 through 18.6. RDQR prints the card input data and diagnostics 
for errors detected. The edited QR input data is written onto the scratch random file 
SCRAND. 


3.3 PRIMARY OVERLAY (L219,2,0)-EOMMOD 
Purpose of EOMMOD 

The primary overlay EOMMOD is called to modify all the matrix equations from the 
L217 (EOM) program according to the instructions specified on the input cards. 

Analytical Steps of EOMMOD 

EOMMOD has the following nine steps: 

1. Read edited input data from SCRAND 

2. Establish pointers for variably dimensioned matrices M], M2, M3, and the active 
control system and sensor data 

3. Check field length available against that required 

4. Establish FET and buffer areas for input EOM file and read the header record if 
the EOM file was generated in the DYLOFLEX system. 


















Table 5.— Routines Called by RDQR 


OVERLAY (L219,1,3) 
PROGRAM RDQR 


RDQR calls 


KEYWRD 
NAMFIL + 
WRITMS * 


+ indicates DYLOFLEX library routine 

* indicates 6600 system library routine 
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If- 

't) 

5 . Call subroutine M 123 to read, modify, and write the Mi, M2, and M3 matrices with 
scalar, replacement, increment, sensor, active control system definition, and body 
axis data 

6. Establish pointers for variably dimensioned matrices M4, M5, and fg and C3 or <f). 

7 . Check field length available against that required 

fsj 

8. Call subroutine M 45 CP to read, modify, and write all the M4, M5, fg , C3, and </> 
matrices with derivative, scalar, replacement, increment, and body axis data 

9 . Print error diagnostics when errors are detected 

Figure 9 is a macro flow chart of EOMMOD. Table 6 displays the subroutines called by 
EOMMOD. 

I/O Devices of EOMMOD 

EOMMOD reads data from the following files: 

EOMTAP Header array plus the matrices Mi, M2, M3, FREQM, M4, M5, 

and C3 (or fg and </>) 

LODTP 2 (Only if sensor equations are being added) the matrices Mi, M2, 

and M3 

SCRAND Edited input EOM data including derivatives, sensor data, and 

active control system data 

EOMMOD writes on two files, the printed output file and EQEOM: 

EQEOM Header array plus the modified matrices Mi, M2, M3, FREQM, 

M4, M5, and C3 (or fg and </>) 

OUTPUT Modified matrices written on EQEOM are also printed. 

3.4 PRIMARY OVERLAY (L219,3,0)-LODMOD 
Purpose of LODMOD 

The primary overlay LODMOD is called to modify all the matrix equations from the 
LOADS program according to the instructions specified on the input cards. 

Analytical Steps of LODMOD 

LODMOD has the following nine steps: 

1 . Read edited input data from SCRAND 
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Figure 9.— Macro Flow Chart of Overlay (L219, 2, 0).—EOMMOD 
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Table 6. —Routines Called by EOMMOD 


OVERLAY 

PROGRAM 


EOMMOD 


(L21 9,2,0) 
EOMMOD 


calls 


FETAD+ 

FETDEL + 

Ml 2 3 calls 


M45CP calls 


ADDSAS 

ADDSEN calls 


BODM12 calls 

INCRM calls 

PRNTM 
READMS * 

REA DTP + 

REPLM calls 

SCALM 

"WRTETP + 

ZEROM 


FETAD + 
FETDEL + 
PRNTM 
READTP + 
SENSOR 
WRITMS * 
IRQL + 


IRQL + 


ADDDER calls 
BODM4 calls 
INCRM calls 
PRNTM 


ADDDIN 
READMS * 
IRQL + 
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Table 6.— (Concluded) 


READMS * 

READTP + 

REPLM calls 

SCALM 

WRTETP + 

ZEROM 


READMS * 
READTP + 
REQFL + 
LOCF * 
MAXO * 


+ indicates DYLOFLEX library routine 
* indicates 6600 system library routine 


IRQL + 
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2. Establish pointers for variably dimensioned matrices Mj, M2, M3, and body axis 
edited data 

3. Check field length available 

4. Establish FET and buffer areas for the input LOADS file and read the header 
record if the LOADS file was generated in the DYLOFLEX system 

5. Call routine M123B to read, modify, and write the Mj, M2, and M3 matrices with 
scalar, replacement, increment, and body axis data 

6. Establish pointers for variably dimensioned matrices M4, M5, C3 and (f>. 

7. Check field length available 

8. Call subroutine M45CPB to read, modify, and write all the M4, M5, C3, and 
matrices with scalar, replacement, increment, and body axis data 

9. Print error diagnostics when errors are detected 
Repeat steps 1 through 8 for each load set. 

Figure 10 is the macro flow chart of LODMOD. Table 7 displays the subroutines called 
by LODMOD. 

I/O Devices of LODMOD 

LODMOD reads data from the following files: 

LODTAP Header array plus the matrices Mj, M2, M 3 , M4, M5 and C 3 or <f>. 

SCRAND Edited loads input data and body axis data 

LODMOD writes on two files: 

EQLOD Modified loads matrices^ including the header array and Mj, M 3 , 

M 3 , M4, M5 and C 3 or $ 

OUTPUT Matrices written onto EQLOD will be printed if requested 

3.5 PRIMARY OVERLAY (L219,4,0)-QRMOD 
Purpose of QRMOD 

The primary overlay QRMOD is called to generate rooting and time solution matrices 
(from the EOM and LOADS equation matrices) for the QR program. 
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Loop on load sets 



Figure lO.-Macro Flow Chart of Overlay (L219, 3, 0)-LODMOD 
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Table 7.— Routines Called by LODMOD 


OVERLAY 

PROGRAM 


LODMOD 


(L219,3,0) 

LODMOD 


calls FETAD + 

FETDEL + 


M123B 

calls 

BOD Ml 2 


calls 

WRITMS * 



I NORM 


calls 

IRQL + 



PRNTM 






READTP 

+ 





REPLM 


calls 

IRQL + 



SCALM 






WRTETP 

+ 





ZEROM 




M45CPB 

calls 

BODM4 


calls 

READMS * 



INCRM 


calls 

IRQL + 



PRNTM 






READTP 

+ 





REPLM 


calls 

IRQL + 



SCALM 






WFTETP 

+ 





ZEROM 
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Table 7.— (Concluded) 


READMS * 

READTP + 

REQFL + 

LOCF * 

+ indicates DYLOFLEX library routine 
* indicates 6600 system library routine 

Table 7 {concluded) 
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Analytical Steps of QRMOD 

QRMOD has the following seven steps: 

1. Read edited input data from SCRAND 

2. Establish FET and buffer area for input EOM file 

3. Read EOM equation matrices and write matrices on SCRAND 

4. If Wagner option was requested, establish matrix pointers, call routine QRWAGN 
for QR’s rooting matrices with Wagner functions, and then go to step 7. 

5. If root option was requested, establish matrix pointers, call routine QRROOT for 
QR’s rooting matrices without Wagner functions, and then go to step 7. 

6. If time history option was requested, establish FET and buffer areas for the input 
LOADS file, read LOADS equation matrices and write matrices on SCRAND, 
establish matrix pointers, and call QRTIME for QR’s time solution matrices 

7. Return to calling program. 

Repeat steps 1 through 6 for each QR set. 

Figure 11 is the macro flow chart of QRMOD. Table 8 displays the subroutines called by 

QRMOD. 


I/O Devices of QRMOD 

QRMOD reads from the following files: 

EOMTAP or EQEOM Mj , M 2 , M 3 , M 4 , M 5 and C 3 or 0 

LODTAP or EQLOD Mj, M 2 , M 3 , M4, M5 and C 3 or 0 

SCRAND Edited input data for QR 

Note: The matrices read from EOMTAP and LODTAP are also temporarily stored on 
SCRAND. 

QRMOD writes on two files: 

QRTAP Matrices which are coefficients of S 4 , S**, S 2 , S 1 and S° 

OUTPUT Matrices written onto QRTAP may be printed by request. 

3.6 DATA BASES 

The programs data bases include input and output files plus internal scratch 
(temporary) storage random file and labeled common blocks. 
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ENTER 



Read edited QR input data from SCRAND 


Establish FET and buffer area for EOM tape 


Read EOM matrices M-, Mo, M^, M^, Mg, and 
Cg or 0 from EOMT AP or EQEOM 


SCRANO 


Write EpM matrices Mj, M2, M2, M^, Mg, and 
C2 or 0 into SCRAND 


l 

SCRAND 


Option 


X Wagner 
Root 

Establish pointers for matrices 
M and S 


SCRAND 


OPTION 


T Wagner 


Cali QRWAGN to 
read EOM matrices 
from SCRAND, generate 
S 4 , S 3 , S 2 , S 1 , S^QR 
matrices, and write 
these matrices on QRTAP 


SCRAND 


Call QR ROOT to 
read EOM matrices 
from SCRAND, generate 
S 2 , S 1 , and S°QR 
matrices, and write 
these matrices on QRTAP 


Establish FET and buffer 
area for LOADS tape 


Read_LOADS matrices 
M2, M^, Wg, and or 0 ^ 
from LODTAP or £OLOD 


W 2 , 


LODTAP 


EQLOD 


Write_LOADS matrices 

M.,. M d , W,-, and or rt I 

onto SCRAND 


Establish pointers for 
matrices M and S 


Call QRTIME to read 
EOM and LOADS matrices 
from SCRAND, generate S 2 
S 1 , S°, and vector QR 
matrices, and write these 
matrices on QRTAP 


SCRAND 


Last 
QR set? 


RETURN 


QRTAP 


QRTAP 


QRTAP 


Figure 11. — Macro Flow Chart of Overlay (L219,4 f O)—QRMOD 
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Table 8.— Routines Called by QRMOD 

OVERLAY (L2 1 9 , 4 , 0) 

PROGRAM QRMOD 

QRMOD calls FETAD + 

FETDEL + 

QRROOT calls PRNTM 

QRMADD 
READMS * 
WRTETP + ■ 
ZEROM 

QRTIME calls PRNTM 

QRMADD 
READMS * 
WRTETP + 
ZEROM 

QRWAGN calls PRNTM 

QRMADD 
READMS * 
WRTETP + 
ZEROM 
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Table 8.— (Concluded) 


READMS * 
READTP + 
REQFL + 
WRITMS * 

+ indicates DYLOFLEX library routine 
* indicates 6600 system library routine 

Table 8 


(concluded) 
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3.6.1 INPUT DATA 


The input data is from two sources, cards and magnetic flies. 

Card Input Data 

For a complete description of all card input formats see section 6.0 in volume I of this 
document (User’s Guide). 


Tape Input Data 

For a complete description of the tape input data see section 6.0 in volume I of this 
document (User’s Guide). 

3.6.2 OUTPUT DATA 

The output data may be of two types, printed and magnetic files. 

Printed Output Data 

For a complete description of the printed output data see section 6.0 in volume I of this 
document (User’s Guide). 


Magnetic Files (Tape or Disk) 

For a complete description of the magnetic file output data see section 6.0 in volume I 
of this document (User’s Guide). 

3.6.3 INTERNAL DATA 

Two methods are used to pass data from one portion of the program to another, labeled 
common blocks and a scratch (temporary) storage random file SCRAND. 

Magnetic File (Scratch Disk File) 

EQMOD uses a random storage scratch disk file for temporary storage of data. All data 
are written on the random disk file SCRAND using subroutine WRITMS. Later, all data 
are read using the subroutine READMS. Table 9 shows the matrices written onto 
SCRAND. The contents of each matrix are then described on the following pages. 


Table 9.— Contents of "SCRAND" 


Matrix Description 


Edited EOM input data 
Edited LOADS input data 


Edited QR input data 


Edited Derivative data 
Edited Sensor data 

Edited active control system definition data 
Anti- symmetric body axis y cQ ^ vector 


EOM matrices 


M. 


M- 


< 


for QR 


m 3 


> 


Mi 


LOADS matrices 


Mi 


M- 


< 


M 3 


for QR time history 


> 




M, 


EOM matrix C 3 for QR time history 
LOADS matrix C 3 for QR time history 


Index 

Name 

Variable 

Length 

EOMA 

LNEOM. 

LODA, 

LODB , . . . 
LODT 

LNLOD (1) 
LNLOD (2) 
LNLOD (20) 

QRA , 

QRB , . . . 
QRT 

LNGR(l) 
LNQR ( 2 ) 
LNQR ( 20 ) 

DERIV 

LNDER 

SENSOR 

LNSEN 

SAS 

LNSAS 

BOD 

LNBOD 

Ml 

LNM ( 1 ) 

M2 

LNM ( 2 ) 

M3 

LNM (3) 

M4 

LNM { 4 ) 

M5 

LNM (5) 

M1B 

LNMB (1) 

M2B 

LNMB ( 2 ) 

M3B 

LNMB ( 3 ) 

M4B 

LNMB ( 4 ) 

M5B 

LNMB (5) 

C3 

LNC3 

C3B 

LNC3B 


36 






Edited EOM Input Data 


SCRAND 
EOMA 

(82+15-NF+NREP+NINC) 
where: NF = Number of frequencies 

NREP = Number of elements to be replaced 
NINC = Number of elements to be incremented data 

Elements: 

Item 1 INEOM, input EOM tape name 

Item 2 INEOMF, file position number where EOM matrices are to be read 

from (after first read, set to zero) 

Item 3 Pointer to header matrix of EOM input tape 

Item 4 Pointer to header matrix of EOM output tape 

Item 5 Pointer to scalar data 

Item 6 Pointer to replacement data 

Item 7 Pointer to increment data 

Item 8 - 37 Header matrix of EOM input tape 

Item 38 - 67 Header matrix of EOM output tape 

Item 68 - (67 + 5 + 3-NF) 

Scalar data 

Item (73 + 3 NF) - (72 + 3- NF + 10- NF + NREP) 

Replacement data 

Item (77 + 9- NF+NREP) - (78 + 9- NF+ 10- NF + NREP + NINC) 

Increment data 

Generation: Program RDEOM 

Edited LOADS Input Data 

SCRAND 

LODA, LODB, .... LODT (up to 20 LOAD sets) 

(83 + 15 • NF + NREP + NINC) 
where: NF = Number of frequencies 

NREP = Number of elements to be replaced 
NINC = Number of elements to be incremented data 


File: 

Index name: 
Dimensions: 


File: 

Index name: 
Dimensions: 
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Elements: 


Item 1 INLOD, input LOADS tape name 

Item 2 INLODF, file position number where LOADS matrices are to be read 

from (after first read, set to zero) 

Item 3 NLDOU, number of output loads 

Item 4 Pointer to header matrix of LOADS input tape 

Item 5 Pointer to header matrix of LOADS output tape 

Item 6 Pointer to scalar data 

Item 7 Pointer to replacement data 

Item 8 Pointer to increment data 

Item 9-38 Header matrix of LOADS input tape 

Item 39 - 68 Header matrix of LOADS output tape 

Item 69 - (68+5 + 3 -NF) 

Scalar data 

Item (74 + 3- NF) - (73 + 3- NF+4 + 6- NF + NREP) 

Replacement data 

Item (78 + 9 -NF + NREP) 

Increment data 


Generation: 


File: 

Index name: 
Dimensions: 
Elements: 
Item 1 
Item 2 
Item 3-32 
Item 33 


Item 34 
Item 35 
Item 36 
Item 37 


Program RDLOD 

Edited QR Data 

SCRAND 

QRA, QRB, QRT (up to 20 QR sets) 

70 x 1 

INEOM, input EOM tape name 

INEOMF, file position number of EOM matrices 

Header matrix of EOM input tape 

ITYPE, indicator for QR process 
= 1, Wagner 

= 2, root 

= 3, time solution 
Wagner function, a 
Wagner function, b 
Wagner function a 
Wagner function, /3 
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Item 38 
Item 39 
Item 40 
Item 41 - 70 
Generation: 


File: 

Index name: 
Dimensions: 


Elements: 
Item 1 
Item 2 


Item 3 


Item 4 


INLOD, input LOADS tape name 
INLODF, file position number of LOADS matrices 
NLDQR, number of loads on LOADS tape 
Header matrix for LOADS input tape 
Program RDQR 

Edited Derivatives Data 


SCRAND 

DERIV 

(21 + 7 -NCS) 

where: NCS = Number of control surface 


NCS, number of control surface 

! x col for symmetric ) 

y col for antisymmetric / 

! z col for symmetric ) 

<f) col for antisymmetric / 


column number of the | y !' freedoms 


column number of the j ^ j ; freedoms 


j d col for symmetric 
I \p col for antisymmetric 


column number of the J ^ j freedoms 


Item 5-(4 + NCS) 5 co i column numbers of the 5 control surface freedoms 

Item (5 + NCS) - (4 + NCS + 6) 

Elements of EOM matrices, six values of 
J x col ) 

\ Ycol / *" irst ^ ree ^ or ^4 and nex f three for Mg 


Item (5 + NCS+6) - (4 + NCS + 12) 

Elements of EOM matrices, six values of 
j z col i 

< > first three for M 4 and next three for M* 

l col J 4 5 


Item (5 + NCS+ 12) - (4 + NCS+18) 

Elements of EOM matrices, six values 



first three for M 4 and next three for M5 



Item (5+NCS+18) - (4+NCS+18+6-NCS) 

For each control surface, NCS, repeat six value of S co j; first three 
for M4 and next three for M5 

Generation: Program RDEOM 

Edited Sensor Data 

File: SCRAND 

Index name: SENSOR 

Dimension: (9 + 2-NM1 + 2-NM2 + 2-NM3) 

where: NM1 = Number of in/out row numbers for Mj matrix 

NM2 = Number of in/out row numbers for M2 matrix 

NM3 = Number of in/out row numbers for M3 matrix 

Elements: 

Item 1 INSEN, input LOADS tape name of sensor 

Item 2 INSENF, file position number where LOADS sensor matrices are to 

be read 

Item 3 NLDSEN, number of loads on sensor matrices 

Item 4-6 Null matrix indicators for three LOADS sensor matrices 

Item 7 NM1, number of in/out row numbers for Mj matrices 

Item 8 NM2 number of in/out row numbers for M2 matrices 

Item 9 NM3, number of in/out row numbers for M3 matrices 

Item 10 - (9 + 2-NM1) 

In row number of sensor matrix M] and out row number of EOM 
matrix Mj 

Item (10+2-NM1)- (9+2- NMl-f 2- NM2) 

In row number of sensor matrix M2 and out row number of EOM 
matrix M2 

Item (10+ 2- NM1 + 2-NM2) - (9 + 2- NMl + 2- NM2 + 2- NM3) 

In row number of sensor matrix M3 and out row number of EOM 
matrix M3 

Generation: Program RDEOM 

Edited Active Control System Definition 

File: SCRAND 

Index name: SAS 

Dimension: (1 + 5 -NS AS) where NSAS = number of Active Control System data 
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Elements: 


Item 1 
Item 2 
Item 3 
Item 4 
Item 5 
Item 6 
Generation: 


NSAS, number of Active Control System data 

ISAS, ith row of augmented matrices 

JSAS, jth column of augmented matrices 

AM1IJ, value of Active Control System element for Mi (I,J) 

AM2IJ, value of Active Control System element for M2 (I,J) 

AM3IJ, value of Active Control System element for M3 (I,J) 

Program RDEOM 

Header Matrix in Edited EOM/LOADS Matrix 


The header matrix contains thirty (30) words 


Word Contents 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 - 20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 . 


7HYDLOFLX 

Program name/version; i.e., L217A1, L219A1, 
Date of run (10H yr/mo/da) 

NDOF, number of degrees of freedom 
NLD, number of load equations 
NPAN, number of panels 
NFREQM, number of frequencies 
of, dynamic pressure 
V, velocity (true air speed) 

(future use) 

Null matrix indicator for 


EOM/LOADS 

Mi, Mi 

M2, M2 
M3, M3 

M 4 , m 4 
m 5 , m 5 

Mg, Mg 
C 2 ,C 2 
c 3> C 3 
f t_ 

l 


Note: If null matrix indicator is zero, the corresponding matrix is null and will not be 

on tape. 

If null matrix indicator is greater than zero, the corresponding matrix does 
appear on tape. 
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Scalar Data in Edited EOM/LOADS Matrix 


Word 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Contents 
Scalar for matrix 


EOM/LOADS 
M 1? M x 
M2, M2 
M 3 , M3 
FREQM 
M4, M4 
M 5 , M 5 

c 2 ,c 2 

C 3) C3 

f 8_ 


Repeat words 5 through 10 for each frequency (NF) times. 
Length is 3 + 5 * NF, where NF = number of frequencies. 


Word 

1 

2 

3 

4 

5 ' 

6 

7 ► 

8 

9 

10. 


Replacement Data in Edited EOM/LOADS Matrix 


Contents 


EOM/LOADS 


Replacement data pointer to/number of 


Word Content 


First 30 bits 

Second 30 bits 

Pointer to 
replacement data 

Number of replacement 
data for matrix 


Repeat words 5 through 10 for each 
frequency (NF). 


Mj, 

M2, M2 
M 3 , M3 
FREOM, 
M4, M4 

m 5 , m 5 

C 2 ,C 2 

c 3> c 3 



(5+6*NFREQM) 


l — Replacement data 


(4+6* NFREQM + NREP) J 

where: NF 

NREP 


= Number of frequencies 
= Number of replacement values 
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Increment Data in Edited EOM/LOADS Matrix 


Word Contents 


EOM/LOADS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


) 


Increment data pointers to/number of 
Word Content 


First 30 bits 

Second 30 bits 

Pointer to increment 
data 

Number of 
increment data 


Repeat words 5 through 10 for each 
frequency (NF) 


Mi.Mj 
M2, M2 
M 3 , M 3 
FREQM, 
M4, M4 

m 5 ,m 5 

c 2 ,c 2 

c 3 >c 3 

f 8_ 


(5 + 6*NFREQM) 


(4 + 6*NFREQM + NINC) J 


Increment data 


where: NF 

NINC 


Number of frequencies 
Number of increment values 


Common Blocks 

Table 10 shows the common blocks used in the program and the overlays in which they 
are defined. 

The labeled common blocks are used for communication between the primary and 
secondary overlays. The block names and contents are described in table 11. 
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Table 1 1.— Common Block Names and Contents 


LABELED COMMON NAME: BODY AX 

DESCRIPTION: Body axis input data. 


VARIABLE 


IBODYA 


DIM. ENG. NOM 


DESCRIPTION 


2 ICOLX 


3 ICOLZ 


4 ICOLT 


5 ICOLY 


6 ICOLP 


ICOLS 


8 BODYA1 


9 BODYVT 


Indicator for body axis 
= 0 , no body axis 
f 0, perform body axis trans- 
formation 

Xcol the column in the matrices 
which are changed by the sym- 
metric body axis transformation 

Zcol the column in the matrices 
which are changed by the sym- 
metric body axis transformation 

®col the column in the matrices 
which are changed by the sym- 
metric body axis transformation 

^col the column in the matrices 
which are changed by the anti- 
symmetric body axis transfor- 
mation 

^ col the column in the matrices 
which are changed by the anti- 
symmetric body axis transfor- 
mation 

^col the column in the matrices 
which are changed by the anti- 
symmetric body axis transfor- 
mation 

“1. angle of attack for body 
axis transformation 

V Tk velocity (true air speed) for 
Dody axis transformation 


NOTE: The "T" heading refers to 
variable type: 

I - Integer 
R - Real 
C - Complex 
L - Logical 
H - Hollerith 
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Table 11.— (Continued) 


LABELED COMMON NAME: HEADER 

DESCRIPTION: Header matrix first nine (9) words. 


NO. 

VARIABLE 

H 

DIM. 


DESCRIPTION 

— 

1 

IWORD1 

H 



The literal DYLOFLEX 

2 

I WORD 2 

H 



Program's name/version which 
the literal "L219A1" 

3 

I WORD 3 

H 



Date of run, YR/MO/DA 

4 

I WORD 4 

I 



Number of degrees of freedom, 
NDOF 

5 

I WORD 5 

I 



Number of load equations, NLD 

6 

I WORD 6 

I 



Number of gradual penetration 
panels, NPAN 

7 

IWORD7 

I 



Number of frequencies, NFREQM 

8 

W0RD8 

R 


q 

Dynamic pressure 

9 

WORD 9 

R 


V T 

Velocity (true air speed) 


is 
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Table 1 1. —(Continued) 











Table 11.— ( Continued) 


LABELED COMMON NAME: NERROR ___ 

DESCRIPTION : Error indicators and accumulators 


NO 


VARIABLE 


DIM. 


ENG. NOMJ 


DESCRIPTION 


NFATAL 


NWARN 


IRR 


Number of fatal errors accumu- 
lated 

Number of warning errors accu- 
mulated 

Error number returned from 
subroutine called 
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Table 1 1 .—(Continued) 


LABELED COMMON NAME: OUTVOL 


DESCRIPTION : 


Output volumes and file position. 



NO. VARIABLE 


NG. NOM 


DESCRIPTION 


IUTEOM 

IFLEOM 


IUTLOD 

IFLLOD 


IUTQR 

IFLQR 



Output EOM tape name, EQEOM 

File position number where EOM 
Matrices will be written 

Output LOADS tape name , EQLOD 

File position number where 
LOADS matrices will be written 

Output QR tape name, QRTAP 

File position number where QR 
matrices will be written 
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Table 1 1.— (Continued) 


LABELED COMMON NAME: PRNTOP 

DESCRIPTION: Print Options 



DESCRIPTION 


Print options for input 
matrices : 

= -999, print all input 
matrices 

= 0 , no input matrix printed 
= I, only matrices of the Ith 
frequency printed 

Print options for output 
matrices : 

= -999, print all output 
matrices 

= 0 , no output matrix printed 
= I, only matrices of Ith 
frequency printed 
= 999, only modified matrices 
printed . 

Checkout print option 
= 0, no checkout print 
/ 0, checkout print. All 

matrices (I/O) and inter- 
mediate results printed. 
Used for debugging 
purposes only. 
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Table 1 1.— (Continued) 



DESCRIPTION 


Number of total degrees of 
freedom (<_ 100) 

Number of degrees of freedom 
read from input tape 

Number of total panels {<_ 50) 

Number of frequencies (<_ 2 0) 

Indicator for type of analysis: 
= 0, symmetric 
^ 0, anti-symmetric 
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Table 1 1.— (Continued) 


LABELED COMMON NAME: SCRANF 


DESCRIPTION:. Scratch random file (SCRAND' 
lengths . 


names and 





DESCRIPTION 


1 NEOM 

2 MAXEOM 

3 LNEOM 

4 NAMEOM 


5 NLOD 

6 MAXLOD 


7 LNLOD 


8 NAMLOD 


9 NQR 

10 MAXQR 

11 LNQR 

12 NAMQR 


13 LNDER 


14 NAMDER 


15 LNSEN 

16 NAMSEN 


17 LNSAS 


18 N AMS AS 


19 LNBOD 


20 NAMBOD 


Number of EOM sets 

Maximum number of EOM sets (<_ 1) 

Length of edited EOM input data 

Matrix name of edited EOM input 
data 

Number of LOADS sets 

Maximum number of LOADS set 

(1 20 ) 

Lengths of edited LOADS input 
data 

Matrix name of edited LOADS 
input data 

Number of QR sets 

Maximum number of QR sets (<^20) 

Length of edited QR input data 

Matrix name of edited QR input 
data 

Length of edited derivative 
data 

Matrix name of edited deriva- 
tive data 

Length of edited sensor data 

Matrix name of edited sensor 
data 

Length of edited active control 
system data definition 

Matrix name of edited active 
control system data definition 

Length of anti-symmetric body 
axis M 3 y col vector 

Matrix name of anti-symmetric 
body axis M 3 y co i vector 
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Table 1 1.— (Continued) 


LABELED COMMON NAME: SCRANF (concluded) 

DESCRIPTION : 

















NO. 

VARIABLE 

T 

DIM. 


DESCRIPTION 

21 

22 

23 

24 

25 

26 

27 

28 

LNM 

NAMM 

LNMB 

NAMMB 

LNC3 

NAMC3 

LNC3B 

NAMC3B 

I 

H 

I 

H 

I 

H 

I 

H 

5 

5 

5 

5 


Length of EOM matrices for QR 

Matrix names of EOM matrices 
for QR 

Length of LOADS matrices for 
QR time history 

Matrix names of LOADS matrices 
for QR time history 

Length of EOM C 3 matrix for QR 
time history 

Matrix name of EOM C 3 matrix 
for QR time history 

Length of LOADS C 3 matrix for 
QR time history 

Matrix name of LOADS matrix 

for QR time history 
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T able 1 1.—( Con tinued) 


LABELED COMMON NAME: TITLE 

DESCRIPTION: Title card 


VARIABLE T 


DESCRIPTION 


1 NTITLE I 

2 ITITLE H (8,4) 


Number of title cards 

Title cards to be printed at 
beginning of printed output 
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Table 11. — ( Concluded) 









BLANK common is used in all secondary overlays and most primary overlays as a 
variable length working area. In general the main program of an overlay calculates the 
area required for arrays in the various subroutines and passes a dimension and first 
word address of each array through the subroutine calling sequence. A description of 
the area used by each overlay is given in figure 12. The BLANK common required 
varies with problem size. Section 6.2, volume I of this document, explains how to 
calculate the core storage required for a particular problem. 
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Field Length per 1000 Octal 


Figure 12.— Overlay Core Requirements and Blank Common Working Area 













4.0 EXTENT OF CHECKOUT 


Five different data cases were assembled to test L219 (EQMOD). The various options 
used are displayed in table 12. The results for each test case were compared against 
hand calculated answers. 
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Table 12.— Options Used in Checkout Data Cases 


Cases 


Options and Maior Paths 

1 

2 

3 

4 

5 

1 

Problem size 

X 

X 

X 

X 

X 

2 

Output tape specification 

X 

X 

X 


X 

3 

Symmetric analysis 

X 

X 

X 


X 

4 

Anti-symmetric analysis 




X 


5 

Body axis transformation 



X 

X 


6 

Equation of motion data 

X 

X 

X 

X 

X 

7 

Sensor 

X 

X 


X 

X 

8 

Derivatives from card 


X 



X 

9 

Derivatives from tape (FLEXSTAB) 

X 



X 


10 

Unsteady derivatives 


X 




1 1 

SAS 

X 




X 

12 

Scale EOM matrices 



X 



13 

Replace EOM matrices 


X 




1 4 

Increment EOM matrices 


X 




15 

Loads equation data 


X 

X 

X 


16 

Scale LOADS matrices 


X 




1 7 

Replace LOADS matrices 



X 



18 

Increment LOADS matrices 



X 



19 

QR data 


X 


X 


20 

Wagner option 


X 




21 

Root option 


X 


X 


22 

Time option 


X 


X 


23 

Diagnostics 
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